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Abstract 

The objective of this study was to analyze the spray characteristics according to the injection duration under ambient pressure 
conditions and the injection timing in the visualization engine. In order to investigate spray behavior, the spray velocity is obtained 
through the PIV method a useful optical diagnostics technology and the vorticity calculated from the spray velocity component. These 
results elucidated the relationship between vorticity and entropy which plays an important role in the diffusion process for early injection 
and the stratification process for late injection cases. In addition, the homogeneous diffusion rate of spray was quantified using the 
entropy analysis based on Boltzmann’s statistical thermodynamics. Using these methods, it was found that the homogeneous mixture 
distribution is more effective as a momentum dissipation of surrounding air than that of the spray concentration to increase the injection 
timing. It is found that the homogenous diffusion rate increased as the injection timing moved to the early intake stroke process, and 
BTDC 60° was the most efficient injection timing for the stratified mixture formation during the compression stroke. 

© 2006 Elsevier Inc. All rights reserved. 
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1. Introduction 

As the environmental problems caused by vehicle 
exhaust emissions have become more severe, exhaust emis¬ 
sion standards and fuel economy regulations have become 
more stringent. For the gasoline engine, the emission of 
C0 2 gas, which is one of the main causes of global warm¬ 
ing, has become a severe problem along with the emission 
of toxic gases such as CO, HC, and NO x . Recently, the 
GDI (gasoline direct injection) engine has been spotlighted 
as the next generation engine that can satisfy the SULEV 
(Super Ultra Low Emission Vehicle) regulations and 
reduce fuel consumption. Thus, studies of a high-pressure, 
vortex-type injector, one of the key parts in the develop¬ 
ment of the GDI engine, have been performed by many 
researchers [19-25]. 
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Zhao et al. visualized the development process of the 
fuel spray of the GDI injector by using a two-dimensional 
Mie scattering method while measuring the SMD (Sauter 
mean diameter) and the velocity of spray to investigate 
the evaporation of fuel [24]. Yamauchi et al. performed a 
numerical analysis on the fuel spray vortex of a high-pres¬ 
sure GDI injector by using a GTT (Generalized Tank and 
Tube) code, and comparing the results with those obtained 
by PDA [25]. In fact, while many other studies have inves¬ 
tigated the characteristics of fuel injector sprays, the exper¬ 
imental difficulties of obtaining detailed information, have 
limited the amount of information available on spray struc¬ 
ture and the characteristics of the GDI injector. 

Yamakawa et al. developed a new fuel spray measure¬ 
ment technology known as, LIF-PIV (Laser induced fluo¬ 
rescence-particle image velocimetry), to measure the 
velocity of fuel spray induced airflow with variations in 
the ambient pressure [8]. They investigated particle diame¬ 
ters, velocity distribution, and the evaporation process of 
the spray. Yuyama [26] and Chikahisa [27] developed a 
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Nomenclature 



S* 

dimensionless entropy value 

An ax 

the maximum intensity in a cell 

M 

number of cells in an image 

h 

image intensity of a cell (0-255) 


integrated value of the image intensity over the 

K 

Boltzmann constant 


whole image area 

CO 

vorticity (1/s) 


new entropy analysis method based on statistical thermo¬ 
dynamics and applied this method to the investigation of 
the propagation process of diesel spray [27,29,30]. How¬ 
ever, these previous studies indicate that there is still much 
uncertainty about the mixing process of gasoline-direct 
injectors. From this point of view, a more exact investiga¬ 
tion of the mixture formation and the correlation between 
spray and induced air is required. 

The present study developed a laser scattering image 
detection system, and measured the macro-scale fuel spray 
characteristics such as spray penetration length and spray 
angle. An entropy analysis method based on the concept 
of statistical thermodynamics has been developed, and 
the mixture formation process was analyzed using this 
method. Cross-correlation DPIV [1-19] has been devel¬ 
oped, and the velocity distribution of the GDI fuel particle 
has been investigated in the visualization engine. And, 
using entropy analysis, the stratification degree near the 
spark plug has been investigated. 

2. Experimental setup and methods 

2.1. Experimental setup and proceedings 

Fig. 1 shows the visualization system used to obtain the 
mixture distribution in the cylinder for simulating intake 
stroke in accordance with the changing intake flow field. 
A Nd:YAG laser (Contiuum Surelite I—10) was used as a 
light source, and the maximum power and wavelength 
was 200 mJ and 532 nm, respectively. The thickness of 
the induced-laser sheet beam was 500 pm. The exposure 
time of the CCD camera, the illumination time of the laser 
beam, and the injection time of the fuel spray were con- 



Optical window 
(cylinder liner) 
(Piston head) 


High pressure chamber 


Double 
Exp. CCD 
Camera 


Encoder signal 
Injector signal 


Image 

proscessor 


Sheet Beam gen. 


DG535 

Injector 

Driver 

Puls* G*n*>atoi 


Half mirror 



Nd:YAG 

Laser 

(Double 

Pulse) 

532nm 


Fig. 1. Schematic diagram of optical visual engine system. 


trolled by a timing board (PC-TIO-10) and the LabView 
program. By controlling the laser trigger signal and delay 
timing, scattering images of fuel spray could be obtained 
through the image board (frame grabber, metio II/digital) 
of a computer by using a high-resolution CCD camera 
manufactured by the Kodak Co. (1008 x 1018, Megaplus 
ES 1.0). Also, we used a pulse generator (DG535, SRS), 
which can be controlled by pico units. 

The particles are illuminated with a 1 mm thin light 
sheet, delivered from an Nd:YAG laser through a combi¬ 
nation of sheet beam generator with spherical and cylindri¬ 
cal lens. Since only a fraction of 1/10 of the maximum 
power of the laser (100 mJ) is employed in the final beam. 
One individual image frame or a sequence of up to 30 
images of the flow can be captured with a CCD camera, 
mounted perpendicular to the light sheet and connected 
to a computer via a frame grabber. The camera requires 
that the laser produces two pulses on every second frame 
in order to achieve optimum image quality. The entire sys¬ 
tem is synchronized using the composite video output sig¬ 
nal extracted from the camera, which represents the input 
of the control circuit. The schematic connections of the sys¬ 
tem components and timing chart are shown in Figs. 1 and 
5, respectively. The specifications of the visualization 
engine and injector are shown in Tables 1 and 2. Table 3 
shows the characteristic specification of the test injector. 


Table 1 

Specification of engine 

Item 

Specification 

Engine 

4 valve DOHC single cylinder 

Fuel 

Gasoline 

Clearance volume 

47.3 cc 

Bore x Stroke 

81 mm x 88 mm 

Compression ratio 

10.3 

Operation speed (rpm) 

600 


Table 2 

Specification of GDI injectors 

Factor 

Injection Spray 


Offset angle 

Swirl 

injector 

pressure pattern 




M 

5 MPa Hollow cone 

0° 

o 

Table 3 





Specification of test injectors 




Type 

0 o (spray angle) 

Injector hole 

Inj. Rate 

Swirl 

CD 

o 

0\ 

O 

0 

0.5 (mm) 

10 (mg) 
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Fig. 2(a) shows the schematic diagram of the injector 
attached angle and optical window. The cylinder liner 
and piston are composed of quartz (fused silica). The fuel 
spray formation process and the mixture distribution in a 
single cylinder visualization engine were observed. A sche¬ 
matic diagram of the laser beam system is shown in Fig. 1 . 
The laser sheet beam was illuminated into the cylinder liner 
or the transparent piston crown. A rotary encoder is 
mounted at the camshaft, and the injection timing and 
the ICCD are synchronized with the encoder signal. The 
spray images were obtained at various crankshaft angles 
during the intake stroke, and then the spray velocity and 
vorticity were calculated using the PIV and vorticity algo¬ 
rithms. Fig. 2(b) shows the schematic diagram of swirl type 
injectors used in this experiment. 

The measurement system of the mixture formation at 
the late injection mode is shown in Fig. 3. To analyze the 
stratification effect near the spark plug at the late injection 
mode, the laser sheet beam was illuminated on a horizontal 
cross-section plane near the spark plug. Optical accesses 
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Fig. 3. Measurement system of the mixture distribution for the late 
injection mode. 
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Fig. 2. Configuration of test engine and specification of the test injector, 
(a) Configuration, (b) swirl injector. 
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were installed at the top of the cylinder head gasket. The 
piston head was modified and transparent. 

The schematic diagram of the visualization region at the 
late injection mode is shown in Fig. 4. The laser sheet beam 
was located 2 mm at the bottom of spark plug. Using this 
system, we analyzed the fuel stratification characteristic 
near the spark plug in accordance with variable injection 
timing as the ignition timing is BTDC 20°. 

2.2. Analysis method 

2.2.1. Cross-correlation particle image velocimetry 

The PIV is now a well established technique for the mea¬ 
surement of instantaneous planar velocity fields and has 
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been reviewed by a number of authors [1-5]. In its simplest 
form, the flow is illuminated with a double pulsed light 
sheet and the positions of tracer particles are recorded with 
a photographic camera viewing normal to the plane of the 
sheet. The mean particle displacement vector in each small 
region of the flow is determined by performing a spatial 
correlation of the particle image field in that region. If 
the image is a simple double-exposure, then autocorrela¬ 
tion analysis is used. If the first and second particle image 
fields can be recorded separately, the more robust cross¬ 
correlation interrogation methods may be used [12,14,18]. 
This gives the important benefits of improved measurement 
dynamic range, better tolerance to velocity gradients 
and seed concentration variations and elimination of 
directional ambiguity. The photographic recording of a 
high-speed time series of moving particle fields for cross¬ 
correlation PIV has been reported for combustion studies 
in constant volume vessels [26] and also in cycle-resolved 
engine studies [32]. It has been shown by practitioners of 
digital PIY that, depending on the details of the flow struc¬ 
ture, reasonable accuracies may be achieved by using inter¬ 
rogation region dimensions as low as 16 x 16 pixels, with 
the condition that particle images occupy no less than 
2x2 pixels [2,4,8,15]. Thus, the possibility exists of gener¬ 
ating velocity measurements on a regular grid of 16 x 16 
locations with a 50% overlap each calculated using an 
FFT-based cross-correlation of a 16 x 16 pixel interroga¬ 
tion region over an image field of 128 x 128 pixels. This 
resolution is adequate for the characterization of the 
large-scale flow structures encountered in typical IC 
engines after inlet valve closure and before the rapid flow 
distortion is encountered in the last stages of compression. 
The available resolution may be used to study a large 
region with a relatively sparse measurement grid or a PIV 
probe system could alternatively be used to measure flow 
to higher spatial resolution or accuracy in smaller regions 
[4,6,13,17]. 

Seeding particles are introduced into the flow in order to 
reveal the motion of the gas. They must be present and accu¬ 
rately follow the flow, while scattering sufficient light to be 
detectable. The inertia of particles affects how well they 
respond to changes in spray and air velocity and in general 
impose the condition that particles are small. This experi¬ 
ments is able to produce (and hence to use in the PIV Sys¬ 
tem) polystyrene particles, whose diameter is about 20 pm. 

The present study recorded particle images exposed by 
each sequential laser beam at separate frames by using a 
CCD camera and calculated the velocity vector using a 
cross-correlation PIV algorithm. The direction of the parti¬ 
cle velocity vector was determined clearly because the 
particle images of each laser pulse were recorded at sepa¬ 
rate frames. Fig. 5 shows the ‘frame-straddling’ method, 
which aligns a laser pulses at each sequential CCD camera 
frames. Fig. 6 shows the cross-correlation PIV algorithm 
using fast Fourier transform (FFT), which is programmed 
to obtain the velocity by calculating the coordinates of the 
location of the maximum cross-correlation parameter. 
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Fig. 5. Schematic diagram of frame straddling and timing chart. 
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Fig. 6. Schematic diagram of the cross-correlation PIV algorithm. 


2.2.2. Verification of PIV algorithm by a standard image 
For the calculation of the velocity field, the PIV uses 
data of area observation instead of that of point observa¬ 
tion, which is used by previous methods. Therefore, the 
PIV can measure the velocity field at a certain instant 
quantitatively, and is widely used in the study of turbulent 
or unsteady flow [1,4]. However, as the experimental appa¬ 
ratus and the analysis algorithm are developed, each 
researcher uses different apparatus and different algorithms 
to increase the accuracy. To establish a standard PIV many 
studies are performed throughout the world. One of the 
studies is the project of JPIV (PIV-STD project) that is held 
by the Visualization Society of Japan. The project evalu¬ 
ated efficiency and accuracy of a PIV system and estab¬ 
lished the PIV standard for analyzing algorithm. 
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The present study evaluated the accuracy of the devel¬ 
oped PIV algorithm with the standard images suggested 
by JPIV. Fig. 7(a) shows one of four sequential reference 
images downloaded by the Visualization Society of Japan, 
and Fig. 7(b) shows the velocity field of the image calculated 
by the program of the present study. The specification of 
standard images is shown in Table 4, and the compared 
results of calculated vector on standard images are shown 
in Table 5. The calculated average velocities show lower val¬ 
ues than the standard with the error of 2%. This is because 
when we divide the image area into small interrogation 
areas, we divide the size of interrogation areas too small. 
Determining the optimal size of the interrogation areas is 
important because the total velocity field shifts down as 
the size of the interrogation area becomes smaller. The pres¬ 
ent study determined the size of the interrogation area by the 




b 


—. ID cm/s 



X-dlrecdonal length scale (mm) 


Fig. 7. Standard image and calculated velocity distribution, (a) standard 
image, (b) velocity distribution. 


Table 4 


Specification of standard images 


Specification 

Value 

Unit 

Image size 

256 x 256 

pixel 

Area 

10 x 10 

mm 

Laser thickness 

2 

mm 

Interval 

0.033 

s 

Average velocity 

8.9 

cm/s 

Particle displacement (average) 

7.5 

pixel/interval 

Particle displacement (standard deviation) 

3.0 

pixel/interval 

Particle displacement (maximum) 

15.0 

pixel/interval 

Number of particle 

4000 

— 

Particle diameter (average) 

5.0 

pixel 

Particle displacement (standard deviation) 

1.4 

pixel 


Table 5 

Compared results of calculated vector on standard images 


Correlating images 

Particle 

displacement 7.5 
(pixel/interval) 

Average velocity 8.9 
(cm/s) 


Value 

Error (%) 

Value 

Error (%) 

1 & 2 

7.402 

1.254 

8.768 

1.473 

2 & 3 

7.352 

1.972 

8.792 

1.212 

3 & 4 

7.384 

1.791 

8.720 

2.014 


usually accepted one-quarter rule, which is suggested by 
Adrian. The rule recommends determining the size of the 
interrogation area as four times the average particle moving 
distance. By the verification using the standard images, we 
can confirm that the velocity field can be obtained by the 
present PIV system with the uncertainty of 2%. 

2.2.3. Principle of vorticity measurement 

The vorticity field, co , is an important Gallilean invariant 
vector field. Each of the vorticity components can be com¬ 
puted from the corresponding measured in-plane velocity 
vector field data using it is the basic definition [6,16]. 

CD = V x u (1) 

where u is the velocity field. This paper is concerned with 
the measurement of one component, denoted by cd z which 
points in the positive Z-direction of an X-Y-Z Cartesian 
coordinate frame. 

dv d u 

cd z = --— 

dX d Y 

where u and v are the velocity components in the V-direc- 
tion and T-direction of the X-Y-Zco ordinate frames, 
respectively. The results derived in this paper apply equally 
well to the remaining vorticity components co x and co y that 
point in the V-direction and T-directions respectively. 

In this study, we used the Ossen vortex algorithm. Using 
the synthetic u and v velocity arrays derived from Ossen 
expression [28], the synthetic vorticity arrays was numeri¬ 
cally generated. Two prevailing techniques were evaluated. 
First, synthetic vorticity was calculated using a simple 
difference equation the standard vorticity algorithm in 
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NASA’s 3-D plotting software, as well as TSI’s post-pro¬ 
cessing software. The data shows (see Fig. 8). 




A L 


Uij -1 1 _|_ Yjlhl YjzhL 


(3) 


where, L means the length of the arrogation area. 

Second, numerical vorticity was calculated using an 
eight-point circulation vorticity equation taken from Reuss 
et al. 


A L 


cd 5 


(2a iy 


U 4 ) 


[u 2 - us) + [u 6 
+0 5 ( — Ug ) + ( U{ ~ Ul ) 

\ — U 2 ) 3r (t /3 — U\) 



Finally, most textbooks on numerical methods point out 
that higher order techniques provide higher accuracy 
because the order of approximation is improved. This does 
not necessarily hold true for noisy data, however, and so 
we have tested on the second order, the central difference 
approximation given by ref [28] in order to ascertain how 
our techniques affected the estimator given here 
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Fig. 8. Schematic of sampling by the vorticity estimator presented in Eqs. 
(3)—(5) [28]. (a) Simple difference, (b) eight-point circulation, (c) second 
order central difference. 


This is one of several second order approaches one 
could adopt. A central difference scheme is the most appro¬ 
priate for comparison to the other first order estimators. 


2.2.4. Entropy analysis 

The entropy analysis of the laser scattering images is 
based on the concept of statistical thermodynamics [31]. 
Boltzmann studied the correlation between the probability 
of particle distribution in the control volume and the 
entropy. If the particle number of N t exists in energy level 
the number of possible combinations ( W) of the total 
particles number N. This is expressed as 


W- N C Nl x • • • x Nm C Nm 

N\ N\ 

- N x \N 2 \N z \---N m \ - UNil 

The entropy might be expressed in the natural logarithm 
of W combination, according to the concept of Boltz¬ 
mann’s statistics. When the Stirling approximation is 
applied to the case of n 1, the natural logarithm of W 
is expressed as 



\nn\ ~ n\nn — n (n 1) 

\nW = \nN\ - ^2 In V 

= N In N -N - N, In N, + ^ N, 
= NlnN - 
S = kln(W) 


= k 


N • ln(N) - Ew • ln (v)i 


(7) 

( 8 ) 


where, S is the entropy and k is a Boltzmann constant. 

If the number of particles in each cell of the image, N h is 
assumed to be proportional to the image intensity of the 
cell, I (ip the entropy can be expressed as 

{ M ^ M 

| • A[ 7 « ■ ln f 7 (')>] ( 9 ) 



where M is the number of cells in the image and a is a figure 
including a Boltzmann constant and the proportional co¬ 
efficient between the number of overall particles and the 
image intensity.If the particle is assumed to be uniformly 
scattered, the average intensity is 



where I t is the integrated value of the image intensity over 
the whole image area. The entropy of a homogeneous 
divergence is then 


Ai = a • A • ln(M) (11) 

The image intensity is divided between 0 and 255, with 
the total image intensity kept constant. So the entropy S 0 , 
obtained from the divided image intensity, is calculated as 

A 0 = a • [I t In (I t ) - PI t • ln(/ MA x)] (12) 
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Where, P is the number of cells that gray light 255 
occupy. 

The normalized entropy S* can be defined as 

C* _ S - So _ h In (/max) ~ Ef {/(■■) • MM1 

Si-S 0 I, • (ln(M) — ln(/,) + ln(/ MA x)} 1 j 

The homogeneity degree of the spray and the diffusion 
phenomena [27,29,30] can be analyzed by normalized 
entropy. 

2.2.5. Verification of entropy analysis 

For the verification of the entropy analysis method, the 
entropy analysis is performed to the images of water parti¬ 
cles coming out of a humidifier. Fig. 9(a) is the image 
upstream near the exit nozzle of the humidifier, and 
Fig. 9(b) is the image downstream 20 cm distance from 
the nozzle. The size of the images is 150x 150 pixels. 
Because the water particles come out from the small nozzle 
of the humidifier and spread to the surrounding area, the 
image downstream shows more homogeneous particle dis¬ 
tribution than the image upstream. From the basic concept 
of the entropy, the entropy at the downstream point should 
be higher than that at the upstream point. The results of 
entropy analysis show higher entropy downstream as 
shown in Table 6. 

Fig. 10 is the comparison of the images of circles with 
different sizes, numbers, sharpness and gray scales but with 
the same total image intensities. The entropy is zero when 
all particles are located at one point, and it increases as the 
particles propagate to the surrounding area. Comparing 
image 1 and image 2 in Fig. 7, image 2 is the image of 
the greater particle propagation, because image 2 has a 



Fig. 9. Comparison of apparently different particle distributions, (a) 
Upstream, (b) downstream. 


Table 6 


Results of the entropy analysis for test images 


Pictures 

Figure No. 

Normalized entropy S* 

Upstream 

Fig. 9(a) 

0.822 

Downstream 

Fig. 9(b) 

0.973 

Image 1 

Fig. 10(a) 

0.288 

Image 2 

Fig. 10(b) 

0.339 

Image 3 

Fig. 10(c) 

0.460 



Fig. 10. Comparison with different size, number, sharpness and gray scale 
of circles in same total image intensity, (a) Image 1, (b) image 2 (c) image 
3 . 


wide spread uniform strength of large circle while image 
1 has circles with particles concentrated at the center of 
the circles. Comparing image 2 and image 3, image 3 is 
the image of the propagation more progressed because it 
has a larger circle. The results of entropy analysis show 
the same results as shown in Table 6. From the entropy 
analysis of Figs. 9 and 10, the normalized entropy calcu¬ 
lated by Eq. (13) has a larger value as the image strength 
is distributed more uniformly and more homogeneously 
throughout the whole area of the image. Therefore, it can 
be verified that the homogeneity degree of the measure¬ 
ment area can be calculated by the entropy analysis, and 
the diffusion strength of the spray can be measured quanti¬ 
tatively from the variation of entropy. 

3. Results and discussion 

3.1. Flow characteristic in the cylinder 

A straight port performance engine block was used in 
this study. Fig. 11 shows the result of the valve lift in accor¬ 
dance with the variable crank angle. 

3.2. Flow velocity distribution in cylinder 

Fig. 12 shows the results obtained for the flow velocity 
in the cylinder in accordance with changing the crank angle 
or valve lift for the intake stroke as the number of rotations 
is 600 (rpm). It calculated the flow velocity from ATDC 90° 
to ATDC 150° by 20 °CA. As shown in Fig. 11, it found 
that the maximum intake flow rate is 130 °CA. From the 
results, the flow entered intake port separates toward both 
walls and in the cylinder occurs at the reverse tumble flow 
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due to the shape of the intake port. The reason for this is 
the performance of the intake port. The velocity of the left 



Crank Angle (°) 

Fig. 11. Variation of the valve lifts in accordance with the crank angle. 


side is trong but that of right side is diminished. In the case 
of the commercial port, that of the right side is strong. 
Accordingly, as the crank angle increases, the velocity of 
the left side is active; then after ATDC 150°, it is inactive. 
But the velocity of the right side diminishes gradually in 
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Fig. 13. Average velocity according to crank angle. 
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Fig. 12. The velocity and vorticity distribution in accordance with crank angle at intake stroke: (a) ATDC 90°, (b) ATDC 110°, (c) ATDC 130°, 
(d) ATDC 150°. 
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accordance with the increasing crank angle. It is observed 
that the reverse collision tumble flow slightly decreased 
and then formed a vortex. As the crank angle increased, 
the vortex forming the interactive collision gradually 


decreased. It expects that the reverse tumble flow due to 
the intake flow has an effect on the homogeneous mixture 
distribution of the fuel droplet accompanied by the flow 
in the cylinder. 
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Fig. 14. Fuel mixture distribution in accordance with crank angle at intake stroke: (a) Injecting timing: 90° ATDC, (b) injecting timing: 110° ATDC, (c) 
injecting timing: 130° ATDC, (d) injecting timing: 150° ATDC. 
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3.3. Velocity in cylinder 

Fig. 13 shows the results of maximum and mean velocity 
in accordance with the increasing crank angle using PIV 
technique as seen in Fig. 6. The strongest part of the veloc¬ 
ity is the middle in the cylinder and then, as the crank angle 
is increased, the mean and maximum velocity in cylinder 
gradually decreases. As shown in Fig. 12, it found that 
the maximum average velocity in cylinder flow is 9 m/s, 
at 110° ATDC; as the crank angle increases, the mean 
velocity gradually decreases. 

3.4. Fuel distribution using Mie-scattering method 

Fig. 14 shows the fuel distribution in the cylinder in 
accordance with changing the injection timing using the 
Mie-scattering method. In the case of early injection 
timing, 90° ATDC, the fuel spay is attached to the wall¬ 
wetting piston crown. If it is injected at 90° ATDC, the 


engine will generate hydrocarbon. Through visual results, 
it found that the injection fitting timing is from 110° ATDC 



Fig. 16. The maximum velocity relative to injection duration. 
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Fig. 15. The spray developing process and spray velocity distribution relative to injection duration ((a)—(c): 2 ms at injection duration, (d)-(f): 4 ms at 
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to 130° ATDC. If it injected the injection timing at 150° 
ATDC, the spray injected from the injector is not involved 
in the reverse tumble flow, as seen in Fig. 12. Therefore, it 
found that if the timing of the injection is retarded after 
ATDC 130°, a homogeneous mixture distribution involv¬ 
ing the intake air does not occur. 


3.5. The velocity distribution of the fuel spray by using 
the PIV technique 

3.5.1. Spray velocity to the free spray 

Fig. 15 shows the velocity distribution of the vertical 
tomogram of the spray at 2 ms and 3 ms from the ignition 
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Fig. 17. The spray velocity and vorticity in accordance with changing injection timing, (a) Injecting timing: 90° ATDC, (b) injecting timing: 110° ATDC, 
(c) injecting timing: 130° ATDC, (d) injecting timing: 150° ATDC. 



































































































































































































































































































590 


K.H. Lee, C.H. Lee / Experimental Thermal and Fluid Science 31 (2007) 579-592 



Time after end of injection (msec) 

Fig. 18. The average vorticity in accordance with changing the injection 
timing. 


of the injection at an ambient pressure of 0.1 MPa. The 
leading edge and central part of the spray shows relatively 
high velocity. Especially, although the velocity scale is 
small, a clear counter-rotating vortex motion is observed 
on the periphery of the spray. In particular, by carefully 
investigating the direction and scale of the velocity vectors 
in the spray periphery, the velocity vectors were found to 
indicate directly the negative axial direction in the axial dis¬ 
tance of 20-30 mm with a small vortex scale in the case of 
the injection duration of 2 ms and 30-50 mm with a large 
vortex scale in the case of the injection duration of 4 ms. 


Fig. 16 shows the maximum velocity in accordance with 
injection duration such as 2 ms, 3 ms, and 4 ms. As the 
injection duration is increased, the maximum velocity after 
the end of the injection in accordance with the elapsing tim¬ 
ing is increased. The maximum velocity to injection dura¬ 
tion is 10 m/s, at 2 ms, 12 m/s, at 3 ms, and 15 m/s, at 
4 ms, respectively. 

3.5.2. Spray velocity and vorticity distribution in the cylinder 

Fig. 17 show the velocity and vorticity distribution of 
spray droplet in accordance with injection timing using 
the cross-correlation PIV technique and second order cen¬ 
tral difference method of the Ossen vorticity algorithm. 
Generally vortex is shown in the left and right but the vor¬ 
tex on right side is not observed in this result due to the 
intake flow. From the Fig. 12, as the velocity is high in case 
of 90-110° ATDC the location of the generated vortex 
moves to the middle direction of in cylinder and disturbs 
the spray velocity distribution. Also it found that the spray 
velocity is involved in the intake flow and the fuel spray is 
dispersed on the in cylinder homogeneously. In particular, 
investigating carefully examine the direction and scale of 
velocity vectors and vorticity strength in the cylinder, it is 
obseved that the velocity vectors indicate directly the neg¬ 
ative axial direction in axial distance of 25-35 mm with 
ATDC 90-130° and about 20 mm with ATDC 150. 

Fig. 18 shows the vorticity of the elapsing timing after 
injection. As seen in Fig. 17, when the injection time is at 
110° ATDC, the vortex size is the largest. And in the case 




Fig. 19. The stratified characteristics relative to injection timing at the compression stroke using the entropy analysis: (a) 40° BTDC, (b) 50° BTDC, (c) 60° 
BTDC, (d) 70° BTDC, (e) 80° BTDC. 
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of 1.6 ms after the end of the injection (AEOI), the vorticity 
strength is highest at 110° ATDC. Finally, we found that 
this condition is optimized when involved in the intake flow 
and spray behavior of the visual engine. 

3.6. The stratification effect of fuel with a compression stroke 
using entropy analysis 

Fig. 19 shows the fuel stratification and mixture proce¬ 
dure using entropy analysis. The study is investigated on 
the entropy analysis of statistical thermodynamics to ana¬ 
lyze quantitatively the diffusion of the mixture distribution 
and homogeneous and heterogeneity degree. The stratifica¬ 
tion effects near the spark plug at the compression stroke 
and set the injection timing from 40° BTDC to 80° BTDC 
at 10° CA intervals were observed. When the ignition tim¬ 
ing was set to 20° BTDC, stratification near the spark plug 
in accordance with the injection timing in the late mode 
was observed and occurred. In the case of the injection 
timing at 50° BTDC and 60° BTDC, the stratification effect 
is higher than in other conditions. Near the spark plug, 
most of the heterogeneity is small scale of the entropy value 
and there are few large clouds. Far away the spark plug, 
the heterogeneities become unclear and the value of that 
is large scale. As the injection timing was slowed and 
speeded up, the stratified effect near the spark plug became 
less. 

4. Conclusion 

A new entropy analysis method based on statistical ther¬ 
modynamics and the cross-correlation PIV measurement 
method has been developed. In this study, we understood 
the optimal conditions and characteristics through analysis 
of the stratification effect and the velocity field measure¬ 
ment of the cylinder. 

By using the Mie-scattering method, entropy analysis, 
and particle image velocimetry, we observed the wall¬ 
wetting of the piston and the effect on the flow field in 
the cylinder in accordance with the fuel distribution and 
mixture homogeneous degree. Through this method, we 
found that the spray velocity is involved in the intake flow 
and the fuel spray is dispersed on the in cylinder homo¬ 
geneously at an injection timing of 110° ATDC. 

Based on the entropy concept of statistic dynamics, a 
method to analyze the homogeneity or heterogeneity 
degree of a field was developed as a tool to investigate phe¬ 
nomena in diffusion process. The homogeneous diffusion 
rate increased as the injection timing moved to the early 
intake stroke process, and BTDC 60° was the most efficient 
injection timing for the stratified mixture formation during 
the compression stroke. From the entropy analysis results, 
it is known that the spray is widely distributed according to 
the elapsed time after injection. Especially, in vaporization 
condition, entropy value is increased in a whole spray area 
from reference. Besides, as ambient pressure increases, the 
fuel distribution becomes heterogeneous condition. How¬ 


ever, the temperature condition in case of the late injection 
is very high. The spray characteristics (entropy value is 
small scale) in case of the injection timing at 50° BTDC 
and 60° BTDC is the state of the heterogeneity which 
appears the low value of entropy. If this experiment is 
achieved in practical operating condition in actual engine, 
the temperature in cylinder become rise. We are expected 
that the mixture distribution near the spark plug change 
the heterogeneity to the homogeneity. 
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